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As batteries become increasingly prevalent in complex systems such as aircraft and elec- 
tric cars, monitoring and predicting battery state of charge and state of health becomes 
critical. In order to accurately predict the remaining battery power to support system 
operations for informed operational decision-making, age-dependent changes in dynamics 
must be accounted for. Using an electrochemistry-based model, we investigate how key 
parameters of the battery change as aging occurs, and develop models to describe aging 
through these key parameters. Using these models, we demonstrate how we can (i) ac- 
curately predict end-of-discharge for aged batteries, and (ii) predict the end-of-life of a 
battery as a function of anticipated usage. The approach is validated through an experi- 
mental set of randomized discharge profiles. 


I. Introduction 


ECENT improvements in battery efficiency and power density, along with lower production cost, have led 

to an increased prevalence in many complex systems, such as smartphones, aircraft, and electric cars. In 
such applications, a critical function is to accurately predict the end-of-discharge (EOD) in order to support 
operational decision-making.t* However, EOD is a direct function of battery age; an older battery will have 
diminished capacity, and, hence, a faster discharge rate.**° Therefore, it is essential to maintain an accurate 
model of battery capabilities over several years of usage, for both EOD and end-of-life (EOL) prediction. 

While some work has been done in understanding the processes and dynamics underlying battery aging,® > 
very little has been done in attempting to estimate battery age, online and in real-time, and using that 
estimate to obtain condition-based predictions of battery discharge and life. Some approaches, such as Saha 
et al.,4 measure the capacity over a sequence of reference discharge cycles (complete discharge at a reference 
current), and fit some empirical expression to the resulting capacity values. With a reference discharge 
cycle, capacity is measured simply by determining how long it takes to reach a voltage threshold that defines 
EOD, and multiplying that by the reference discharge current to obtain a capacity rating in Ah. Capacity 
measured in this way is dependent on the reference current, because, at a given age, a battery will have 
a lower Ah rating when the reference current is increased, and higher one when the reference current is 
decreased, i.e., the relationship between discharge current and discharge time is nonlinear. Further, if the 
battery is allowed to rest sufficiently in between discharges, the apparent capacity will increase.' Thus, 
there are two significant problems with such an approach. First, measuring capacity in this way is not an 
accurate description of battery age, since the Ah rating will depend on the current (and other factors such as 
temperature), unless such reference conditions can be kept constant (which in an operational environment 
is not practical). Second, reference discharges are not available in real battery usage; real usage is complex 
and often unpredictable. 
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In contrast, our goal is to develop an approach to estimate battery age that can be used in everyday online 
usage, and does not require reference cycles. Our previous work towards that end used an unscented Kalman 
filter (UKF) to estimate jointly the battery state and aging parameters online® using an electrochemistry 
model.' However, the estimation problem in this case is difficult because minor changes in aging parameters 
over the course of a single cycle can easily be hidden in the data, and corrected simply by updating the state 
estimates and keeping the parameter estimates static. Essentially, the dynamics of the aging parameters are 
at a much slower time scale than that of the battery states, so it is difficult to estimate them jointly with the 
states. Instead, in this work, we separate the estimation task into two phases, a real-time phase in which we 
estimate only the battery states, assuming aging parameters are known, and a cycle-time phase, in which the 
aging parameter estimates are updated after every cycle. We assume here that the usage context is that of 
a battery that undergoes distinct charge and discharge cycles, with variable loads, where aging parameters 
can be updated based on all the data in a given discharge cycle. This kind of usage context can be found, 
for example, in electric aircraft, where the aircraft battery is charged pre-flight, and discharges a significant 
amount during flight before landing.? 1° 

The paper is organized as follows. Section II presents the model-based prognostics framework and the 
prognostics architecture. Section III reviews the electrochemistry-based modeling for battery discharge. 
Section IV presents the battery aging models. Section V describes the estimation approaches. Section VI 
discusses end-of-discharge and end-of-life prediction. Section VII concludes the paper. 


II. Model-based Prognostics 


We formulate the EOD and EOL prediction problems using a model-based prognostics framework. In this 
section, we formulate the prognostics problem, according to.!! We then provide a computational architecture 
for battery EOD and EOL prognostics. 


A. Problem Formulation 
We assume the system model may be generally defined as 


x(k + 1) = £(k, x(k), 0(k), u(k), v(k)), (1) 
y(k) = h(k, x(k), A(k), u(k), n(k)), (2) 


where & is the discrete time variable, x(k) € R” is the state vector, 0(k) € R”® is the unknown parameter 
vector, u(k) € R™ is the input vector, v(k) € R”” is the process noise vector, f is the state equation, 
y(k) € R™ is the output vector, n(k) € R™ is the measurement noise vector, and h is the output equation.* 
The unknown parameter vector 0(k) is used to capture explicit model parameters whose values are unknown 
and time-varying stochastically. 

The goal of prognostics is to predict the occurrence of some event FE that is defined with respect to the 
states, parameters, and inputs of the system. We define the event as the earliest instant that some event 
threshold Tg : R"* x R™ x R™ — B, where B £ {0,1}, changes from the value 0 to 1. That is, the time of 
the event kg at some time of prediction kp is defined as 


ke(kp) &inf{k @N: k > kp A Te(x(k), 0(k), u(k)) = 1}. (3) 


























In this paper, we have two separate prognostics problems, and two corresponding system models, one for 
EOD and one for EOL. For EOD prediction, the event EF represents EOD and is determined by a voltage 
threshold; the battery is considered to be at EOD when the voltage is less than the low voltage limit. In 
this case, we will assume there are no unknown parameters (i.e., 9 = @). For EOL prediction, the event F 
represents EOL and is determined by a capacity threshold for a reference discharge; the battery is considered 
to be at EOL when the capacity is less than the given lower capacity limit. 


B. Prognostics Architecture 


In the model-based prognostics architecture,'?:!3 there are two sequential problems, (7) the estimation prob- 
lem, which determines a joint state-parameter estimate p(x(k), O(k)[Y RP ) based on the history of observa- 
tions up to time k, denoted as Xie and (7) the prediction problem, which determines at prediction time kp, 





*Bold typeface denotes vectors, and nq denotes the length of a vector a. 
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Figure 1. Battery EOD and EOL prognostics architecture. 
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using p(x(k),O(k)[YR"), p(Urp), P(Vep), and p(@x,), the probability distribution p(ky(kp)|¥;”). Here, 
U;,, denotes the future system inputs from kp on, V;,, denotes the future process noise values from kp on, 
and ©; denotes the future unknown parameter values from kp on. 

In this paper, we will assume @;, is known exactly, Vx, is zero, and U;, is known exactly, in order 
to focus on validation of the aging models and the performance of the associated predictions. Dropping 
these assumptions will be addressed in future work. A general framework for dealing with all sources of 
uncertainty in prognostics is discussed by Sankararaman et al.!4 

The overall combined EOD/EOL prognostics architecture is shown in Fig. 1. In discrete time k, the 
system is provided with inputs u(k) (current) and provides measured outputs y(k) (voltage). There are two 
models used for the two different prediction problems: one for EOD (with subscript EOD), and one for EOL 
(with subscript EOL). The age parameter estimation block estimates the states for the EOL model, xgox(k), 
based on the data from the previous discharge cycle. In the context of the EOD model, xgoz,(k) are the 
aging parameters and assumed to be constant for a given discharge. The state estimation block estimates 
states for the EOD model, xgop(k), using the latest estimates of the aging parameters. EOD prediction 
computes the probability distribution of the EOD time, kgop, using the latest state estimate and aging 
parameter estimates. The aging rate parameters, Ogoi(k) for the EOL model are estimated based on the 
estimated aging parameters over past discharge cycles, and parameterize how quickly the aging parameters 
change in time. EOL prediction computes the probability distribution of the EOL time, kgop, using the 
current estimates of the aging and aging rate parameters. 

We next describe the EOD and EOL models, followed by descriptions of the estimation and prediction 
algorithms for these models. 








III. Battery Discharge Modeling 


A battery is a collection of electrochemical cells that convert between chemical and electrical energy. Each 
cell consists of a positive electrode and a negative electrode with electrolyte. In this paper, we focus on Li-ion 
cells. The electrolyte enables lithium ions (Lit) to diffuse between the positive and negative electrodes. The 
lithium ions insert or deinsert from the active material depending upon the electrode and whether the active 
process is charging or discharging.!® 

In this work, we use an electrochemistry-based model developed for Li-ion cells in our previous work.! 
Unlike other electrochemistry-based models that rely on complex partial differential equations that are un- 
suitable for online estimation and prediction algorithms, this model uses only ordinary differential equations, 
and is fast enough for real-time use. In this section, we briefly summarize this model and describe its key 
features. 

The voltage terms of the battery are expressed as functions of the amount of charge in the electrodes. 
Each electrode, positive (subscript p) and negative (subscript n), is split into two volumes, a surface layer 
(subscript s) and a bulk layer (subscript b). The differential equations for the battery describe how charge 
moves through these volumes. The charge (q) variables are described using 


ds,p = tapp + Gbs,p (4) 
qb,p = —4s,p + tapp — tapp (5) 
db,n = —Gbs,n + tapp — tapp (6) 
ds,n = —tapp + Gbs,n; (7) 


where igp, is the applied electric current The qq, term describes diffusion from the bulk to surface layer for 
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electrode i, where i = n or i = p: 


Ere eet (8) 


dbs,i = D 


where D is the diffusion constant. The c terms are Li-ion concentrations: 








On = qb, (9) 
Ub,i 

C= (10) 
Usyi 


Here, cy; is the concentration of charge in electrode 7, and v,,; is the total volume of charge storage capability. 
We define v; = vp; + Us, Note now that the following relations hold: 





dp = As,p T Yb,p (11) 
Gn = Is,n as db,n (12) 
qo = s,p T Yb,p + ds,n + db,n- (13) 


We can also express mole fractions (2) based on the gq variables: 


Gi 








i= qmax’ (14) 
qs,i 

tsi = gman (15) 
Syt 
bd 

m= Oe, (16) 
i 


where g™** = q, + qd, refers to the total amount of available Li-ions. It follows that x, + x, = 1. For Li-ion 


batteries, when fully charged, x, = 0.4 and x, = 0.6. When fully discharged, x, = 1 and x, = 0.'° 

The overall battery voltage V(t) consists of several electrochemical potentials. At the positive current 
collector is the equilibrium potential Vy,,. This voltage is then reduced by V;,,, due to the solid-phase ohmic 
resistance, and V,,,, the surface overpotential. The electrolyte ohmic resistance then causes another drop 
V.. At the negative electrode, there is a drop V,,, due to the surface overpotential, and a drop V,,,, due to 
the solid-phase resistance. The voltage drops again due to the equilibrium potential at the negative current 
collector Vy,,. These voltages are described by the following set of equations: 








RT 1-25; 
Vu, = Up + — In ( “s ) + Vint; (17) 
nF sii 
dss 2x,k(1 — 2;) 
4s, yeti See a 
Vint i — nF oa Aj, k (con 1) (2x; r= =) ’ (18) 
Vo — tappRo, (19) 
RT . Jj 
Ve= Fo tresinh (3) ; (20) 


a 
5, 
Jig = ki(1 — 5,4)" (25,4) *; 

V = Vu» — Von — Ve — Van — Va 


1 

( 

,P n,n? ( 

V3 = (Vo —Ve)/To ( 
( 

( 


J, 
22 
23 
24 
25 
26 


) 
) 
) 
) 
ie = (Vip > Vid ete ) 
Van = (Vin — Vin) / Tan: ) 


Here, Uo is a reference potential, R is the universal gas constant, T is the electrode temperature (in K), n is 
the number of electrons transferred in the reaction (n = 1 for Li-ion), F' is Faraday’s constant, J; is the current 
density, Jio is the exchange current density, and k; is a lumped parameter of several constants including a 
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Figure 2. Battery voltage during 1 A discharge cycle. 


rate coefficient, electrolyte concentration, and maximum ion concentration. Vint; is the activity correction 
term (0 in the ideal condition). We use the Redlich-Kister expansion with N, = 12 and N,, = 0 (see'). The 
7 parameters are empirical time constants (used since the voltages do not change instantaneously). 

The state vector, input vector, and output vector of the EOD model are defined follows: 


T 
xnop(t) = [ae db,p dbsn Ys,n V3 Vid Vee ’ (27) 


u(t) = [éap| (28) 
y(t) = [Vv]. (29) 


Parameter values for a typical Li-ion cell are given in Daigle et al.1 

An example discharge cycle is shown in Fig. 2. The cell starts fully charged at 4.2 V and is dishcarged 
at 1 A. Around 8000 s, the voltage hits the lower voltage limit, 3.2 V, which defines EOD. At this point the 
load is removed and the cell voltage recovers back to the equilibrium potential. Here, we use 3.2 V as the 
voltage threshold defining EOD. In this case, the capacity is computed as the discharge time (7630 s) times 
the applied current (1 A), yielding in this case 2.12 Ah. 

Note that capacity can only be measured consistently for a discharge cycle at reference conditions, since 
measured capacity is a nonlinear function of the load and environmental conditions. In this work, reference 
conditions are defined by a 1 A discharge at room temperature. 


IV. Battery Aging Modeling 


Generally speaking, battery aging results in a loss of capacity and an increase in internal resistance.” 
This is due to several physical processes such as solid electrolyte interface layer growth, which contributes to 
the resistance increase, diffusion stress, which leads to changes in diffusion properties and contributes to the 
capacity loss, and irreversible parasitic side reactions, which result in a loss of active Li-ions and contribute 
to the capacity loss. These aging phemonena have mainly been studied under constant loading conditions, 
but here, we want to estimate aging parameters based on real-life usage. 

Battery behavior can also change over time due to faults, but here, we limit to normal degradation 
and aging processes. The effect of aging on a normal discharge cycle is shown in Fig. 3. Clearly, capacity 
decreases, since the time of EOD gets pushed earlier and earlier. Further, resistance is increased, accounting 
for the drop in voltage during discharge and also to the earlier EOD time. The diffusion parameter is also 
increasing, resulting in a slower diffusion rate, and thus the apparent state-of-charge (that observed at the 
electrode surface) decreases faster than with a normal diffusion rate, also contributing to the decrease in 
EOD time. Equilibrium potential after discharge is higher because less charge is drawn from the battery, 
due to the earlier EOD time and the decrease in diffusion rate (causing the 3.2 V limit to be reached sooner). 

For the EOL model, we focus on how aging parameters change in time as a function of usage. Following 
our previous work,° we select as the aging parameters g™**, capturing the decrease in active Li-ions, and 
Ro, capturing the resistance increase, but, unlike our previous work, also consider the diffusion constant D, 
to capture changes in diffusion properties that can affect capacity. In the context of the EOD model, these 
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Figure 3. 
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Figure 4. Battery voltage during 1 A discharge cycle halfway to EOL. 
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Figure 5. Battery voltage during 1 A discharge cycle at EOL. 


parameters are constant, but in the context of the EOL model they are states: 
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(30) 


Figs. 4 and 5 show how well changes to these parameters can capture the aging process, shown here 
for halfway to EOL and at EOL (where EOD, and thus measured capacity, is half of nominal). Although 
there are some differences, the time of EOD and the equilibrium potential after discharge are captured very 


These variables change as a function of battery usage, i.e., as a function of the applied current, depps” UE 
We assume that the following models descibing their dynamics: 


(31) 
(32) 
(33) 


where the w parameters are the aging rate parameters, which are constants. If w,, wr, wp, and future 
values of tapp are known, then we can predict future values of xgoL, and given that, determine whether the 
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Figure 6. Estimates of aging parameters over discharge cycles. 


capacity as measured at reference conditions would be below the limit. Here, we use a 50% of the nominal 
capacity as the threshold defining EOL. 

Our experiments consist of a sequence of charge and discharge cycles at randomized load. Since capacity 
can only be measured for a discharge at reference conditions, reference discharge cycles have been performed 
at regular intervals in order to obtain a ground truth EOL. Note that in practical usage scenarios, refer- 
ence discharges will not be available and capacity can only be estimated from the estimates of the aging 
parameters, XEOL- 


V. Estimation 


In this section we describe the estimation procedures. As shown in Fig. 1, there are three estimation 
problems: (7) state estimation for the EOD model, (ii) state estimate for the EOL model, and (7) parameter 
estimation for the EOL model. We describe each of these in turn. 


A. State Estimation for the EOD Model 


State estimation for the EOD model has been developed in previous work,!’!! so here we only summarize 
the approach. Here, we have a standard state estimation problem, with the model described by that in 
Section III and the states given by xgop. Since the model is nonlinear, a nonlinear filter is needed. For this 
purpose, we use the unscented Kalman filter (UKF),'*!® due to its high accuracy and low computational 
cost. The UKF uses the concept of sigma points, which are deterministically selected samples that are 
selected to represent the first two statistical moments of a distribution. The unscented transform is the 
procedure that selects sigma points from a given distribution. By passing sigma points through a nonlinear 
function, an approximation to the first two statistical moments of the original distribution passed through 
that function is computed. The filter first updates the current state estimate using the state equation, then 
corrects the estimate based on the observations, through the use of the sigma points. Details of the algorithm 
and its application to prognostics are provided elsewhere.!> 2° 


B. State Estimation for the EOL Model 


For the EOL model, the state estimation problem is to estimate xpoz, i.e., g™°*, Ro, and D. In this case, 
we assume these parameters are approximately constant over a given discharge cycle, and so use system 
identification techniques to estimate their values for a given discharge cycle. For this, we use the Nedler- 
Mead algorithm for multi-dimensional unconstrained nonlinear minimization. For given values of the aging 
parameters, the model is simulated and the sum of the mean-squared error is computed between the predicted 
voltage and the observed voltage. The optimization algorithm searches over values of the aging parameters 
to minimize this error. The algorithm is very fast, and finds a solution in under 20 s on average. For optimal 
results, the discharge cycle data must include the rest period after reaching EOD, in order to correctly 
characterize q™**. 

This procedure is repeated after every discharge cycle. As the battery ages, the estimated values of these 
parameters change in time. Fig. 6 shows the results of this procedure on a battery undergoing random walk 
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Figure 8. Estimates of aging rate parameters over discharge cycles. 


discharge cycles. The estimates are fairly noisy, however a clear trend is observable: g™** decreases with 
age, R, increases with age, and D increases with age. 

For given values of the aging parameters, we can simulate a discharge cycle at reference conditions to 
compute a corresponding value of estimated capacity. Fig. 7 shows the estimated capacity as a function of 


the estimated aging parameter values after each discharge cycle. Capacity is most sensitive to changes in 
qh, 


C. Parameter Estimation for the EOL Model 


For estimation of the aging rate parameters w,, wr, and wp, we try to fit values of these parameters to the 
estimated values of the aging parameters. In a discrete-time form, the equation describing g™**(k) is 


Ghe“(k + 1) = g™*(k) + walk) |dapp() |dt(k). (34) 


Note that the discrete time points are after every discharge cycle, and, hence, are not evenly spaced, i.e., 
the time interval dt(k) is not constant. Between discharge cycles, the applied current is variable, so for 
liapp(k)| we take an average of the current for the kth discharge cycle (and the preceding charge cycle). The 
formulation is equivalent for R,(k) and D(k). 

Again, we use the Nedler-Mead algorithm, using the estimated values of the aging parameters from the 
initial value of k, ko, up to the current value. After each new estimate of the aging parameters, we update the 
aging rate parameter estimates based on all available aging parameter estimates up to the current discharge 
cycle. With more data, the estimate becomes more accurate and converges. The estimates of the aging rate 
parameters are shown in Fig 8. The parameter estimates converge after about 100-200 cycles. 

Fig. 9 shows the predicted values of the aging parameters based on the aging models using final estimated 
values of the aging rate parameters, against their estimated values. These aging models describe the evolution 
of the aging parameters as the battery ages fairly well. 
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Figure 9. Estimates and predictions of aging parameters over discharge cycles. 
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Figure 10. Battery voltage during random walk discharge cycle halfway to EOL. 


VI. Prediction 


There are two different prediction problems to consider, EOD prediction and EOL prediction. We use 
the same algorithm in both cases, which is basically to simulate the given model forward for given intial 
state and parameter values, and hypothesized future inputs and process noise. The core procedure is shown 
as Algorithm 1. If probability distributions are given, then this algorithm can be used on realizations of the 
prediction inputs within an algorithm such as Monte Carlo sampling. 

For both EOD and EOL prediction, we assume that the future usage is known, and process noise is neg- 
ligible. The case where future usage is described by a probability distribution is covered in other works.!!!8 
For EOD prediction, we assume there are no unknown parameters (the values used for the aging parameters 
are the most recently estimated by the EOL model state estimation algorithm), and for EOL prediction, we 
assume the aging rate parameters are constant and their values given by the most recently estimated values 
by the EOL parameter estimation algorithm. 

As seen in Figs. 2-5, when the current to be applied is known, accurate predictions of EOD can be 
computed even at t = 0. An open loop prediction for a random walk discharge is shown in Fig. 10, and is 
also very accurate in predicting EOD in the open loop; EOD is predicted to be only 5 s later than it actually 
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Figure 11. EOL predictions. 


occurs. Predictions at later times are only more accurate. However, as the battery ages, there is a bias that 
appears in the EOD predictions, and EOD is overestimated. As a result, capacity estimates will be biased 
and capacity will be overestimated. This means that additional parameters need to be considered to change 
as a function of aging, in order to more accurately characterize the EOD time at higher battery ages, since 
the three selected aging parameters are not sufficient to eliminate that bias. 

For EOL prediction, at any prediction point, we have the current estimated values of the aging parameters 
and aging rate parameters. We can use these and the EOL model to predict the future values of the aging 
parameters and to see when these future values of the aging parameters will result in the capacity being less 
than the limit. Results are shown in Fig. 11. The true EOL is determined to be at the 844th discharge cycle. 
At each prediction point, we use as the initial state the estimated aging parameters at that cycle and use 
the estimated aging rate parameters computed at that cycle. Since the EOD prediction is biased at higher 
ages, the capacity estimates are also biased and the EOL predictions are biased. With a better parameter 
fit at higher ages, this bias can be eliminated. 


VII. Conclusions 


In this paper we proposed a new framework for EOD and EOL prediction for Li-ion batteries, proposing 
new aging models. Procedures for estimation and prediction are described. The overall approach can estimate 
the current battery state, perform EOD prediction using up-to-date aging parameters, and perform EOL 
prediction using aging models and estimated aging rate parameters. 

Here, we selected a critical subset of battery parameters and observed how they change in age through 
system identification procedures applied over any given discharge cycle. Based on the EOL prediction results, 
we found that additional parameters must be considered as aging parameters in order to better characterize 
EOD at higher ages. This would result in better (unbiased) capacity estimates and correspondingly better 
EOL predictions. In addition, no filtering was performed on the aging parameter estimates, but this would 
result in more precise capacity estimates and EOL predictions. Uncertainty in the states, parameters, and 
future usage was also not considered here, and is also a subject of future work. 
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